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Abstract 

Contrairement a une opinion largement repandue, le probleme des 
trois corps possede une solution analytique. Cette solution fut decouverte 
en 1909 par Sundman. Nous presentons dans cet article les idees de 
base et l'histoire de cette solution. 

Entgegen einer weitverbreiteten Meinung ist das Dreikorperproblem 
analytisch losbar. Diese Losung wurde 1909 von Sundman gefunden. 
Die ihr zugrundeliegenden Ideen und ihre Geschichte werden in ein- 
facher Form dargestellt. 
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1 Le probleme des trois corps 



L'etude du probleme des trois corps a une longue histoire [Barrow-Green 
1997, Diacu et Holmes 1996]. Comme beaucoup de problemes qui ont suscite 
un grand interet, il se formule tres facilement mais sa simplicite apparente 
cache une phenomenologie tres riche. Considerons un ensemble de n par- 
ticules ponctuelles qui exercent des forces gravitationnelles l'un sur l'autre. 
Admettons qu'elles ont les masses rn.j, ou i = 1, 2, . . . , n et que leurs positions 
sont decrites a l'aide des vecteurs r*j(t) en fonction du temps t. Les equations 
newtoniennes du mouvement s'ecrivent 

d 2 ri(t) " miTiij fj(t)-fi(t) 
n^i — = — G > 9 ■ , _; ; ; : — Kr ; * = l,2,...,ra (1) 

ou G est la constante gravitationnelle. Le probleme des trois corps s'obtient 
comme cas special n = 3 du probleme general des n corps ([!]). En admet- 
tant des conditions initiales pour les positions r i)0 := n(0) et les vitesses 
Vifl := dfi/dt(0) on s'interesse a calculer les positions fi(t), a partir des 
equations differentielles (|1|), pour tous les temps t. Ainsi pose, le probleme 
est mathematiquement completement defini et sa solution n'utilise que des 
techniques mathematiques relevant d'un cours de mecanique classique.0 

Quelques solutions particulieres sont bien connues. Euler (1765) traite 
le cas ou les trois corps sont alignes. Lagrange (1772) resout le cas oil 
les trois corps (dont un de masse negligeable) sont aux sommets d'un tri- 
angle equilateral rigide. Finalement, dans la solution de Moore (1993) et 
de Chenciner et Montgomery (2001) trois corps de masse egale parcourent 
un huit. Pour des details complementaires, on peut consulter [Montgomery 
2001]. 

De part a sa difficulte technique considerable, le probleme des trois corps 
suscite des commentaires souvent approximatifs et parfois erronnes. Par ex- 
emple, on lit souvent que ce probleme est insoluble et on attribue a Poincare 
la decouverte d'un tel theoreme. Cette insolubilite, reliee a la divergence 
des series que Ton utilisait pour resoudre les equations (ffl), aurait stimule 



3 Des exemples physiques classiques du probleme des n corps sont evidemment fournis 
par le mouvement des planetes autour du soleil, ou par le passage d'une comete dans le 
systeme solaire. Des exemples plus recentes considerent les orbites des satellites artificials, 
ou des planetes recemment decouvertes autour d'autres etoiles que le soleil. 
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Poincare pour developper ses nouvelles methodes de nature topologique en 
mecanique analytique [Abraham-Marsden 1978, xvii]. Sous cette forme sim- 
ple, ces on-dits ne sont que partiellement vrais (et done, en toute rigueur, 
entierement faux). En effet, Sundman [1909] a donne, dans un sens tres 
precis, une solution exacte et generate du probleme des trois corps. Elle 
s'ecrit a l'aide de series convergentes pour toute valeur (reelle) du temps t. 

On peut se demander pourquoi un tel resultat n'est connu qu'aux specialistes. 
Cette note, destinee a une audience non specialised dans les methodes mod- 
ernes de la mecanique analytique, presentera dans la section 2 un tres bref 
resume historique dans lequel nous tenterons de rectifier les erreurs de cer- 
tains commentaires et d'en expliquer l'origine. Dans la section 3 nous presenterons 
un expose des idees essentielles de la solution de Sundman, en suivant large- 
ment un article pedagogique de Saari [1990]. Des conclusions dans la sec- 
tion 4 termineront cette note. En particulier, nous souligenerons le fait que 
la solution de Sundman, malgre son interet theorique fondamental, n'est 
pas capable de fournir des renseignements pratiques sur le comportement 
physique du systeme. Quelques annexes d'intention pedagogique fournissent 
des complements mathematiques d'une nature un peu plus technique. 

2 Quelques remarques historiques 

L'histoire du probleme des trois corps, a la fin du 19e siecle et au debut 
du 20e siecle, est etroitement liee a l'histoire du journal Acta Mathematica, 
fonde et redige par G. Mittag-Leffier a Stockholm. Pour plus de details sur 
les relations entre Mittag-Leffier et Poincare, on peut consulter [Nabonnand 
1999a, 1999b]; pour une vision d'ensemble de l'histoire du probleme des trois 
corps a la fin du 19e siecle, on peut consulter [Barrow-Green 1997] et [Diacu 
et Holmes 1996]. 

Mittag-Leffier a fait ses etudes en Allemagne, chez K. WeierstraB a Berlin, 
chez E. Schering a Gottingen, et en France, chez C. Hermite a Paris. A 
cette epoque, caracterisee par des nationalismes gratuits et des tensions 
politiques franco-allemandes, les contacts et la connaissance reciproque des 
mathematiciens frangais et allemands etaient difficiles. Mittag-Leffier, re- 
connu en particulier pour ses resultats sur des fonctions meromorphes qui 
portent aujourd'hui son nom [Knopp 1971], fut nomme en 1881 professeur 
a Stockholm [Nabonnand 1999b, 71]. Vue la situation politique qui rendait 
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delicate la publication des travaux des mathematiciens frangais dans les jour- 
naux allemands (et reciproquement), il tente de tirer profit de ses bonnes 
relations avec les mathematiciens frangais et allemands en creant un nou- 
veau journal, nomme Acta Mathematica, qui publiera principalement des 
articles en frangais et en allemand, dans le but avoue d'offrir un forum in- 
ternational oil des travaux d'importance peuvent etre publies et lus dans les 
deux pays sans heurter les sensibilites politiques. En mars 1882, il demande 
a Poincare de "faire le succes"Q du journal et de lui envoyer pour publica- 
tion ses memoires sur les fonctions fuchsiennes dont des enonces de result at s 
furent deja publies dans les Comptes Rendus. Poincare accepte et son travail 
[Poincare 1882] inaugure le premier volume des Acta en 1882. II sera suivi 
d'une serie de neuf autres ceuvres dans les 10 premiers volumes des Acta. 
L'entreprise de Mittag-Leffler est soutenue par ses collegues scientifiques et 
le journal acquiert rapidement une bonne reputation internationale, [Nabon- 
nand 1999b, 92]. Cependant, la situation financiere du journal reste delicate, 
et les Acta Mathematica ne peuvent demarrer que grace aux contributions 
de la fortune personnelle de Mittag-Leffler et des subventions du parlement 
suedois (en 1895, la reduction de la subvention par le parlement de Stockholm 
jettera les Acta Mathematica dans une grave crise financiered]) . 

Pour toutes ces raisons, Mittag-Leffler cherche a augmenter le prestige 
de son journal des son lancement. Une telle opportunity se presente a 
l'occasion du 60e anniversaire du roi Oscar II de Suede qui, ayant etudie 
les mathematiques lui-meme, prend un interet personnel dans les progres des 
mathematiques, et accepte de financer un prix pour la solution d'un probleme 
mathematique important. L'histoire et les evenements autour du prix du roi 
Oscar II ont ete decrits en detail ailleurs [Andersson 1994, Barrow-Green 
1994, Barrow-Green 1997, Diacu 1996, Diacu et Holmes 1996, Nabonnand 
1999a] et nous nous concentrons ici sur les aspects en lien direct avec la so- 
lution du probleme par Sundman [1909, 1913] une trentaine d'annees plus 
tard. 

Le prix est annonce, en frangais et en allemand, dans le volume 7 des 
Acta (et egalement dans de nombreux journaux internationaux [Nabonnand 
1999b]) et quatre questions sont proposees. Void la premiere: 

" 1. Etant donne un systeme . . . de points materiels qui s'attirent 

4 Lettre de Mittag-Leffler a Poincare du 29 mars 1882 
5 Lettre de Mittag-Leffler a Poincare du 20 avril 1895. 
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mutuellement suivant la loi de Newton, on propose, sous la suppo- 
sition qu'un choc de deux points n' ait jamais lieu, de representer 
les coordonnees de chaque point sous forme de series procedant 
suivant quelques fonctions connues du temps et qui convergent 
uniformement pour toute valeur reelle de la variable. 

. . . la solution etendra considerablement nos connaissances par 
rapport au systeme du monde . . . Lejeune-Dirichlet a communique 
peu de temps avant sa mort a un ge.ome.tre . . . qu'il avait decouvert 
une methode de I'integration des equations differentielles . . . [et] il 
etait parvenu a demontrer d'une maniere absolument rigoureuse 
la stabilite de notre systeme planetaire. Malheureusement nous ne 
connaissons rien sur cette methode . . . On pent pourtant supposer 
que cette methode etait basee . . . sur le developpement d'une idee 
fondamentale et simple ..." [Mittag-Leffler 1885] 

La suite est connue: Poincare envoie[] le 17 mai 1888 son memoire celebre 
[Poincare 1890] ,[] inscrit sous la devise Nunquam praescriptos transibut sidera 
fines, destine au concours pour le prix du roi Oscar II. Apres deliberation 
et de multiples voyages de Mittag-Leffler, le comite forme de Mittag-Leffler, 
Hermite et WeierstraB [Nabonnand 1999b, 178-181] decide de lui octroyer 
le prix, qui consiste en une medaille d'or et la somme de 2500 couronnes.f] 
Selon le reglement, le memoire gagnant doit etre publie dans les Acta Math- 
ematica. En preparant cette publication, Poincare fournit sur la demande de 
Mittag-Leffler des notes explicatives sur son memoire (ce qui le rallonge d'une 
centaine de pages). En juillet 1889, suite a une question de l'editeur associe 
des Acta, E. Phragmen, Poincare decouvre une erreur importante dans son 
memoire.0 Une grande partie des resultats que Poincare a cru etablir, entre 
autres sur la stabilite du systeme planetaire, ne sont plus valables. Travail- 
lant sur les corrections necessaires, il decouvre ce qu'on appelle aujourd'hui 

6 Tout laisse croire que les quatre questions avaient ete redigees de maniere a pouvoir 
interesser Poincare. En effet, Mittag-Leffler lui rappelle en juillet 1887 que le delai final 
pour envoyer un manuscrit pour le prix est le l er juin 1888 et ajoute: "Si vous veuillez 
envoyer quelque chose c'est guerre probable que quelqu'un vous depassera." 

7 Cet article ne traite que du probleme restreint des trois corps, c.a.d. les corps sont 
admis a ne se deplacer que dans un plan et m 3 <g; rai,m,2. 

8 Pour comparaison, le salaire annuel de Mittag-Leffler a Stockholm en 1881 etait de 
7000 couronnes. En 1894, 1 couronne etait equivalente a 1,40 francs frangais. 

9 Lettres de Poincare a Mittag-Leffler du 16 juillet et du ler decembre 1889. 
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les points homoclines et met en evidence le premier exemple d'un systeme au 
comportement chaotique. Au moment de la decouverte de l'erreur, le vol- 
ume contenant le memoire errone est imprime et quelques exemplaires sont 
deja sortis de l'imprimerie. Afin d'eviter un scandale, Mittag-Leffler decide^ 
de rappeler le volume complet, de le detruire et de l'imprimer a nouveau. 
Poincare doit assurer le cout de cette seconde impression, d'un montant de 
3585 couronnes et 63 0re. Finalement, son memoire du prix du roi Oscar 
II lui servira de base pour son grand ceuvre sur les nouvelles methodes de 
la mecanique celeste [Andersson 1994, Barrow-Green 1994, 1997, Diacu et 
Holmes 1996, Nabonnand 1999b]. 

Ajoutons ici quelques commentaires: 

1. On lit souvent dans les articles de vulgarisation que le probleme des 
trois corps n'est pas resoluble. Formellement, un systeme algebriquement 
integrableQse caracterise par l'existence de certaines quantites Ik = Ik{ri(t), . . . ,f n (t)), 
appelees integrates premieres, qui restent constantes tout au long de 
revolution du systeme. Par consequent, elles ne dependent que des valeurs 
initiales des positions r it0 et des vitesses Vi t o, i = 1, . . . , n, c'est a dire 

h = h (n,o, • • • , r n ,o, #1,0, • • • , #n,o) (2) 

Si Ton connait autant d'integrales premieres Ik que de variables independantes 
(on en a 6n pour le probleme des n corps) et si les Ik sont des fonctions suff- 
isamment simples (c.a.d. algebriques) des r^o et t^o, le probleme de resoudre 
les equations differentielles ([]]) se reduit a trouver la solution d'un systeme 
d'equations algebriques. Ceci est un probleme considerablement plus simple 
que celui de resoudre les equations ([I]) directement. 

On connait 10 integrates premieres independantes pour le probleme des 
n corps. Elles peuvent s'ecrire comme : 

(i) les coordonnees R du centre de masses du systeme, 

(ii) les composantes de la quantite de mouvement totale P, 

(iii) les composantes du moment cinetique du systeme L et 

(iv) l'energie totale E. 

Les equations (jl]) sont des equations differentielles du second ordre, on a 
done 2 • 3n = 6n variables independantes. On peut toujours eliminer 12 = 
10 + 2 variables du probleme, grace aux 10 integrales premieres et 2 par 

10 Lettre de Mittag-Leffler a Poincare du 5 decembre 1889. 

n La terminologie technique definie dans cet article sera mise en gras. 
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deux symetries decouvertes par Jacobi en 1843.0 En somme, on reste avec 
un systeme d'equations differentielles a 6n — 12 = 6(n — 2) variables. II est 
clair que le probleme des n = 2 corps est integrable mais il faudrait trouver 
d'autres integrales premieres pour rendre les cas n > 2 integrables. Or, ceci 
est impossible: Bruns a demontre en 1887 que toute integrale premiere qui 
est une fonction algebrique des positions r*j(t) et des vitesses dr*j(t)/dt est 
une fonction des dix integrales deja connues [Bruns 1887]. Poincare [1890] a 
generalise ce theoreme aux fonctions uniformes (voir annexe A) des positions 
et des vitesses.0 La non-existence d'autres integrales premieres etablit la 
non-integrabilite du probleme des n > 3 corps. Pourtant, ceci n'implique pas 
qu'une solution exprimable a l'aide des series n'existait pas. Le theoreme de 
Bruns et Poincare montre seulement que certaines methodes algebriques sont 
insuffisantes pour resoudre le probleme. 

2. II est bien connu qu'il existe des series perturbatives capables de 
representer les solutions du probleme avec une tres grande precision. Pour 
simplifier considerons le cas n = 3. Imaginons qu'un des trois corps est 
le soleil et les deux autres des planetes. En premiere approximation, on 
neglige les forces gravitationnelles entre les deux planetes, parce que leur 
masse est beaucoup plus petite que celle du soleil.0 Dans cette approxima- 
tion, l'eq. (|l|) se decompose en deux problemes a deux corps et s'integre. 
Ensuite, on rajoute les forces entre les planetes en les traitant comme une 
petite perturbation de la solution approximative obtenue auparavant. On 
trouve une meilleure approximation "proche" de la premiere. En repetant 
cette procedure, on obtient des expressions pour les r*j(t) sous forme d'une 
serie. Le cas des mouvements periodiques presente des difficultes techniques 
qui furent resolues a la fin du 19e siecle par l'utilisation des series de Lindst- 

12 Une variable s'elimine en considerant une des positions ou des vitesses comme une 
variable independante (au lieu du temps t) et l'autre s'elimine a l'aide de Y "' elimination 
des nceuds" qu'on trouvera dans les livres sur la mecanique celeste. 

13 Painleve a generalise ceci en 1897/1900 en demontrant qu'il n'existe aucune integrale 
Jfe algebriquement independante qui soit une fonction algebrique/uniforme des vitesses. 
Aujourd'hui, des theoremes comme celui de Ziglin [1983] fournissent des conditions 
necessaires pour l'integrabilite. On peut en deduire que le probleme des trois corps, au 
voisinage de la solution speciale de Lagrange (1772), ne possede pas de systeme complet 
d'integrales premieres qui soient des fonctions meromorphes des positions et des vitesses 
[Tsygintsev 2000]. 

14 Par exemple, la masse du soleil est environ mille fois celle de la planete la plus grande, 
Jupiter. 
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edt et Gylden.0 Poincare [1890] demontre la divergence generique des series 
de Lindstedt, mais souligne aussi que 

". . . les considerations qui precedent n'enlevent rien au merite 
pratique des developpements de M. Lindstedt. lis ne convergent 
pas; done Us ne peuvent donner une approximation indefinie; 
mail Us peuvent donner assez rapidement une approximation tres 
grande et tres suffisante pour les besoins de la pratique. 

Je serais desole d 'avoir jete quelque discredit sur ces series 
. . . parce que je regarde la methode de M. Lindstedt comme I'une 
des meilleures qui soit connues. " [Nabonnand 1999b, 193 note 4] 

3. Le developpement recherche pour le prix du roi Oscar II n'est pas 
une serie du type perturbatif. On y cherche un developpement convergent 
des coordonnees f£(t) sous la forme Y^Lo a vfif) u - H semble connu depuis 
longtemps que les solutions devraient etre de cette forme, comme il apparait 
dans une lettre de Weierstrafi a Mittag-Leffler de 1883 

"Angenommen nun, es gehe die Bewegung der Art vor sich, 
dafi niemals zwei Punkte zusammentreffen, so sind [fi,i = 1, . . . , n] 
eindeutige analytische Funktionen vont, nicht blofifiir reelle Werte 
dieser Grofie, sondern auch fur alle komplexen, in denen die 
zweite Koordinate ( der Faktor von i ) dem absoluten Betrage nach 
unter einer gewissen Grenze liegt. . . . Man wird schwerlich a pri- 
ori die Bedingungen ermitteln konnen, die erfiillt sein miissen, 
damit niemals zwei Punkte zusammentreffen konnen, man wird 
vielmehr dieselben als erfiillt vorraussetzen miissen . . . Poincare 
hat . . . [dieses] Theorem ebenfalls hergeleitet, wenigstens unter 
Vorraussetzung des Newtonschen Gesetzes und daraus die Fol- 
gerung gezogen, es sei moglich, die Koordinaten aller Punkte in 
konvergierende Reihen der Form Y^Lq a uf{t) u zu entwickeln, wo 
ip(t) eine bestimme Funktion von t ist. Dies ist leicht einzusehen. 
. . . Aber man erhdlt auf diese Weise nicht Aufschlufi daruber, ob 
die gemachte Vorraussetzung erfiillt ist oder nicht und es ist auch 
die Form, in der sich die Ausdriicke der Koordinaten darstellen, 

15 Un exemple simplifie de la methode de Lindstedt sera expose en annexe C. 
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nicht der Natur der zu beschreibenden Bewegungen angemessen. " 
[Nabonnand 1999b, 119]0 

Apres avoir regu de Mittag-Leffler une copie de cette lettre, Poincare remar- 
que qu'en effet la forme <p(t) = tanh(at/2) n'est utilisable que dans le cas 
ou il n'y a pas de collision, mais dans ses propres etudes il travaille avec la 
fonction ip(t) = tanh(as(t)/2), ou s(t) est une fonction monotone de t encore 
a determiner. Selon Poincare :[[] 

"Les solutions de ce probleme [de determination de s(t)] sont 
en nombre infini . . .11 est clair que dans chaque cas particulier, 
il faut choisir la plus zweckmaBig.0 Or je ne crois pas que dans 
le cas de la Mecanique Celeste celle que j'ai donnee soit la plus 
zweckmafiig, je crois qu 'il y a mieux a trouver. " 

II apparait que la possibility des collisions entre les corps est l'obstacle prin- 
cipal pour d'etablir une solution analytique du probleme. 

4. En 1909, Sundman [1909] demontre qu'une solution sous forme de 
series convergentes pour tout temps t existe pour le probleme des trois corps. 
Voici l'enonce de son theoreme: 

"Si les constantes des aires dans le mouvement des trois corps 
par rapport a leur centre commun de gravite ne sont pas toutes 
nulles, on peut trouver une variable r telle que les coordonnees des 
corps, leurs distances mutuelles et le temps soient developpables 
en series convergentes suivant les puissances de r, qui representent 
le mouvement pour toutes les valeurs reelles du temps, et cela 
quels que soient les chocs qui se produisent entre les corps. " [Sund- 
man 1909, 3] 

Sur l'invitation de Mittag-Leffler, Sundman [1913] publie une exposition 
d'ensemble de ses travaux dans les Acta Mathematica. Constatons qu'aucune 
trace de ce resultat ne se trouve dans la correspondance entre Mittag-Leffler 
et Poincare. Le travail de Sundman [1907,1909,1913] fournit la solution 

16 L'orthographe de l'allemand est mis a jour. 
17 Lettre de Poincare a Mittag-Leffler du 22 mai 1883. 
18 utile, pratique; en allemand dans roriginal. 
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recherchee par la question 1, posee par Weierstrafi pour le prix du roi Oscar 
II une trentaine d'annees auparavant, au moins pour le cas n = 3.0 

Nous verrons dans la section prochaine dans quelle mesure les idees de 
Weierstrafi et de Poincare sont proches de celles de Sundman. 

3 La theorie de Sundman 

Nous presentons les idees principales de la solution de Sundman du probleme 
des trois corpsQ Nous allons nous servir du probleme des deux corps afin 
d'illustrer les idees et nous mentionnons les generalisations necessaires pour 
le probleme des trois corps.Q 

II est bien connu que dans le probleme des deux corps, le mouvement est 
dans le plan perpendiculaire au vecteur du moment cinetique^] L. Comme 
L est conserve, son orientation et done aussi le plan du mouvement des deux 
corps sont fixes. II est pratique d'utiliser des variables complexes z = x + iy 
afin de decrire le mouvement dans le plan (voir annexe A pour la terminologie 
des nombres complexes). 

Afin d'en donner un exemple simple, considerons le mouvement d'une 
planete autour du soleil. Dans la figure [l|, nous montrons les orbites de Mars 
et de Jupiter, vues depuis la terre. Leur mouvement apparent semble tres 
complique et il pourrait sember difficile de decrire et comprendre d'un point 
de vue physique une telle orbite. Bien entendu, nous savons aujourd'hui que 
cette difficulte n'est que apparente et qu'en realite, les planetes sont en orbite 
autour du soleil. Ces orbites sont en bonne approximation des cercles, dont 
les rayons et les periodes sont rassembles dans le tableau: 

19 Sundman a trouve la solution pour le cas n = 3. Presque un siecle plus tard, et a 
l'aide de techniques differentes, Wang [1991] a pu obtenir un resultat analogue pour toute 
valeur de n. 

20 Cette section est inspiree de Particle pedagogique de Saari [1990]. 

21 Pour une presentation mathematiquement complete, dans une notation moderne, voir 
[Abraham et Marsden 1978, Siegel et Moser 1971]. 

22 Pour une particule a la position r, de vitesse v et de masse m, le vecteur du moment 
cinetique s'obtient comme L = mvAf, ou A est le produit vectoriel. Pour un systeme de n 
particules, chacune avec un moment cinetique Li, oui= 1,2, ... ,n, le moment cinetique 
total du systeme est L = L\ + L 2 + . . ■ + L n . 
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Figure 1: Mouvement de Mars (orbite interieure) et de Jupiter (orbite 
exterieure) vus depuis la terre (au centre). 
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A l'aide des variables complexes, nous pouvons exprimer les orbites circulaires 
des planetes, vues depuis le soleil, sous la forme 



ZT (t) = 2e 2wit ; z M (t) = Se wit ; zj(t) = 10e^ /6 



(3) 



Le mouvement relatif a la terre s'obtient par un simple changement de repere. 
Par exemple, pour l'orbite apparente de Mars, on a 



z(t) = z M (t) - z T (t) 
= 3 e nit - 2 e 27Tit 
= 2 + e^(3-4cos(7rf)) 



(4) 



et geometriquement, l'orbite est un limagon. De meme, pour Jupiter, on 
obtient un limagon relatif a un cercle tournant 



z(t) = zj(t)-z T {t) 

= 2 e" 57ri '/ 3 + e^ 6 (lO - 4 



llTT 

cos [ 1 



(5) 
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Une telle description geometrique est deja utilisee par Ptolemee et les as- 
tronomes de son ecole. Dans notre terminologie, il admet que les orbites 
apparentes sont de la forme 

parent (t) =aie fel ^ + a 2 e" 2 ^ (6) 

ou les ai, a2, £>i, 62 sont des constantes a determiner. Le premier terme decrit 
le mouvement d'un planete selon le "deferent" et le second decrit le mou- 
vement additionnel selon "l'epicycle". Plus systematiquement, le travail 
des astronomes avant Copernic (Kopernikus) revient a chercher une orbite 
planetaire sous la forme 

00 

^apparent (^) ^ ] Qj & J (7) 
3=1 

et de trouver les aj et bj, c.a.d. d'effectuer une analyse de Fourier d'une 
orbite planetaire. II est clair qu'en incluant suffisamment de termes, on peut 
arrive a un tres bon accord avec les observations. Bien que ce jugement soit 
completement recursif et anhistorique, on ne peut s'empecher de comparer 
la simplicity des expressions lorsque l'on choisit comme origine du repere 
le soleil (|3p et la complication, entierement artificielle, d'expressions comme 
celles qui apparaissent dans @. Un "bon" choix de variables permet de 
simplifier considerablement l'analyse des orbites apparemment compliquees 
de la figure |l|. 

De meme, pour le probleme des trois corps, c'est un nouveau choix de 
variables qui va permettre a Sundman de contourner un obstacle majeur 
a la comprehension de ce probleme. La theorie de Sundman utilise le fait 
(voir annexe A) que la valeur du rayon de convergence d'une serie est reliee 
aux singularites des fonctions complexes. Or, dans le cadre du probleme des 
trois corps, Interpretation de telles singularites est particulierement simple : 
En effet Painleve a demontre en 1895 que pour n = 3 toute singularity des 
fonctions fi(t) par rapport a t correspond soit (i) une collision entre deux 
corps soit (ii) une collision entre les trois corps.0 Afin d'etablir l'existence 
d'une solution sous forme d'une serie convergente, Sundman [1907] part de 
ce theoreme et progresse en deux etapes: 

23 Pour n > 5, des singularites qui ne proviennent pas de collisions simples peuvent 
apparaitre. Pour n = 4, on ne sait pas s'il existe d'autres singularites que des collisions 
ou non. 
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1. Eliminer les collisions binaires entre deux particules. 



2. Exclure les collisions ternaires entre toutes les trois particules. 

Dans le cas des collisions binaires, on peut se restreindre au probleme a deux 
corps.0 En effet, dans ce cas l'orbite a la forme d'une ellipse qui s'ecrit en 
coordonnees polaires sous la forme r{6) = a /{I — ecos9) ~a(l + ecos#) ou a 
est le demi-grand axe et e est l'excentricite de l'orbite. Dans le cas r = 
et lorsque l'excentricite e tend vers 1, la variable angulaire 6 passe d'une 
valeur proche de a une valeur proche de 2ir. Pour obtenir une solution 
analytique, cette singularity en 6 doit etre eliminee. Pour cela, on decrit le 
mouvement dans le plan a l'aide d'une variable complexe z = z(t) = x + iy. 
L 'equation du mouvement s'ecrit alors sous la forme 

d 2 z z , 

^=-- 3 >r=\z\ (8) 

On obtient une regularisation du probleme par le changement de variables 
[Sundman 1909, eq. (11)] 

w = w(s) = z{tfl 2 , ds = ^ (9) 

Avec ces nouvelles variables, l'equation du mouvement s'ecrit 

d 2 w(s) h 

~ds^ ~ 2 W{S) = ° (10) 

ou la constante h correspond a l'energie totale.0 

La solution de l'equation ( PH| ) est pour h ^ 0: w(s) = w cosh s \ + 

Wi sinh (J h/2s^J ou U7 0j i sont des constantes. Elle est sans singularite et se 
developpe en serie convergente en s pour tout |s| < oo. 

Pour le probleme des trois corps l'idee reste la meme mais le temps t est 
remplace par une nouvelle variable to [Sundman 1909, 1913] 

dt = (l - e- r ^) (l - e~ r ^) (l - e"^) du (11) 



24 Lors d'une collision binaire, le 3e corps est loin des deux autres qui se rencontrent. 
Au voisinage de la collision, son influence peut done etre negligee. 

25 Les calculs detailles sont exposes dans l'annexe B. Nous montrerons que ce changement 
de variables est tres naturel d'un point de vue physique. 
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oil 7*0,1,2 sont les distances entre les trois corps et I est une constante. 

Le cas des collisions ternaires semble etre fournie par le fait qu'une col- 

— * 

lision entre les trois particules n'est possible que si le moment cinetique L 
s'annule.0 En effet, la condition L ^ est facile a controler a partir des 
conditions initiales0 et done elimine toute singularite en rj(t) pour toute 
valeur reelle de t. Cependant, un examen plus precis revele une complication 
technique dont nous allons discuter brievement, avant de revenir a la solution 
de Sundman. 

Pour que le rayon de convergence du developpement en serie du rayon 
vecteur fi(t) so it non nul, il est necessaire qu'il n'y ait pas de singularite 
au voisinage de l'axe reel. Or ceci n'est pas trivial a etablir meme dans le 
cas n = 2. Dans ce cas, si l'excentricite e < 1, l'equation de l'orbite s'ecrit 
(equation de Kepler) 

r(v) = a(l — ecosv) , t — v — esin-u (12) 

Bien que cette solution soit reguliere pour les temps reels, la fonction r{t) a 
des singularites complexes aux points 

Vk = 2irk + iarcosh (1/e) (13) 

(oil k est un entier arbitraire) qui sont proches de l'axe reel pour e suffisam- 
ment proche de 1. Mais comme la valeur de e depend de celle du moment 
cinetique L du probleme des deux corps, on peut montrer lorsque L ^ 0, 

r > r min > (14) 

oil la valeur de r min depend de L. 

Sundman [1907] a pu demontrer un resultat analogue pour le probleme 
des trois corps, c.a.d. que si L ^ 0, on a 

\r i (t)-r j (t)\>C(L)>0 (15) 

26 Cette proposition semble avoir deja ete etabli par Weierstrafi, mais Sundman [1907] 
fut le premier a publier la demonstration, voir aussi [Siegel et Moser 1971]. 

27 Si une collision ternaire se realise, on a L = 0. Dans ce cas et au voisinage d'une 
collision, Sundman trouve "que les corps se meuvent de telle maniere . . . [qu'ils] ou bien 
forment de plus en plus un triangle equilateral ou bien se rangent de plus en plus en ligne 
droite." [Sundman 1907] Ce type de configuration s'appelle configuration centrale et est 
aussi present dans le cadre du probleme a n corps. 
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Figure 2: Transformation conforme de la variable complexe u (a gauche) vers 
le plan complexe r = tanh (ttuj/4) (a droite). 

oil C(L) est une constante. Ayant ainsi etabli que pour L ^ 0, il existe une 
borne inferieure pour les distances entre les corps, Sundman considere les 
temps formellement complexes et montre que le developpement en serie des 
rayons vecteurs ne presente aucune singularite dans un ruban |Im u\ < B(L) 
oil B(L) est une constante connue dependant du moment cinetique L et u est 
defini par l'eq. (|Tl|). If suffit alors d'utiliser les transformations conformes 
considerees par WeierstraB et Poincare 

t = tanh(vrcj/4 J B) (16) 

qui projette le ruban |Im u\ < B sur l'interieur du cercle \r\ < 1 (voir figure |] 
pour le cas B = 1). Dans la partie gauche, un morceau du plan complexe en 
io est presente. En exprimant les rayons vecteurs = ^(r) en fonctions de 
£(t), il n'existe aucune singularite pour \r\ < 1. Comme r = ±1 correspond 
k t = ±oo, on a une representation des positions sous forme d'une serie 
uniformement convergente pour toute valeur \t\ < oo. Le probleme des trois 
corps a done trouve sa solution sous la forme initialement conjecturee par 
WeierstraB. 

4 Remarques finales 

La solution de Sundman ne fait appel qu'aux techniques d'analyse standard 
des fonctions complexes, deja connues a son epoque. Dans ce sens, elle 
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confirme l'attente de WeierstraB, formulee lors de l'annonce du prix du roi 
Oscar II, que la solution devrait etre basee sur des idees simples. On peut 
voir aussi que par rapport a WeierstraB et Poincare, Sundman ajoute deux 
nouvelles idees: (i) la regularisation des collisions binaires a l'aide de son 
changement de variables flTTD ce qui contourne la difficulte de caracteriser 
le cas des collisions binaires et (ii) sa demonstration de separation minimale 
des trois corps lorsque L ^ 0. Pourtant, la forme meme de cette solution 
ne fournit aucune information sur le comportement du systeme. Au con- 
traire, le theoreme de Sundman met en evidence que la propriete d'analyticite 
ne peut pas servir a caracteriser l'une ou l'autre propriete qualitative d'un 
systeme avec trois corps. En particulier, aucun renseignement sur la stabilite 
a des temps tres longs n'est fourni. De plus, ces series sont inutiles pour des 
calculs numeriques. Certes, elles convergent, mais le taux de convergence 
est tellement faible qu'un calcul purement numerique, utilisant des tech- 
niques de resolution directe des equations de mouvement (||) sans faire appel 
aux developpements des solutions en series, est plus fiableQ WeierstraB et 
Poincare, ont-ils pressenti ce comportement en constatant que la variable r 
etait "nicht der Natur der zu beschreibenden Bewegung angemessen"^\ et de 
ne pas etre "zweckmdftig" ? 

Initialement, et surtout apres la publication, de son article dans VActa 
[Sundman 1913], le travail de Sundman fut recu avec grand interet,Q mais 
ce travail parait avoir ete tres vite oublie [Barrow-Green 1997]. II n'est pas 
entierement evident si cet oubli soit du au succes enorme des methodes qual- 
itatives de Poincare pour l'etude des sytemes dynamiques ou que quelques 
unes des idees contenues dans les travaux de Sundman puissent encore au- 
jourd'hui stimuler de nouvelles directions de recherche -P 1 ] Siegel et Moser 
[1971] consacrent tout un chapitre de leur traite de la mecanique celeste 
aux travaux de Sundman [1907,1909,1913] et soulignent leur importance. En 

28 Ceci vient du fait que le premier changement de variable (fill), de la forme u> ~ t 1 / 3 , 
ralentit la dynamique. Enfin, le second changement illustre en figure |^ concentre toute la 
dynamique au voisinage de r ~ ±1. 

29 pas conforme a la nature du mouvement a decrire 

30 En 1913, Sundman recevait le prix Pontecoulant de l'Academie des Sciences, dont la 
valeur a ete doublee a cette occasion [Barrow-Green 1997]. 

31 Dans le contexte de la mecanique quantique, rintegrabilitc du cas n = 2 mene a 
une symetrie dynamique, que Pauli en 1927 a utilise dans son traitement dc l'atomc 
d'hydrogene. Existe-t-il une trace de la solution de Sundman en mecanique quantique, en 
tenant compte de la non-integrabilite du cas n = 3 ? 
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revanche, Abraham et Marsden [1978] ne mentionnent qu'en passant les "oeu- 
vres classiques de Sundman (1913)". 

Notons encore que la solution de Sundman est constructive. En relation 
avec le debat sur V intuitionisme de Brouwer [Diacu 1996], il n'est pas sans 
interet de remarquer que meme une approche entierement constructive ne 
donne pas toujours des resultats directement utiles. Cet exempleQ est une il- 
lustration excellente de la difficulte a definir des notions mathematiques d'une 
telle maniere que des resultats profonds peuvent etre obtenus. Dans le cadre 
de l'histoire du probleme des trois corps, dont nous devons la comprehension 
au progres initie par Poincare [1890], et les debuts du journal Acta Mathemat- 
ical cela fournit egalement une illustration parfaite des occasions excellentes 
pour des decouvertes majeures qu'une approche scientifique reellement inter- 
nationale, depassant entre autres des contraintes politiques, peut offrir. 
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Annexe A. Sur les fonctions complexes 

Nous rappelons quelques faits elementaires sur les fonctions analytiques f{z) 
d'une variable complexe z [Knopp 1976]. Un nombre complexe z = x + iy 
est caracterise a l'aide de deux nombres reels x, y. On peut le representer 
comme etant un point (x, y) dans un plan, voir figure |3]a. z redevient reel si 
y = 0. On appelle x = Re z la partie reelle de z et y = Im z la partie 
imaginaire de z. L'unite imaginaire i satisfait a la condition i 2 = —1. Le 
nombre z = x — iy s'appelle le conjugue complexe de z. Alternativement, 

32 La solution satisfait a toutes les demandes de rigueur mathematique, selon les stan- 
dards etablies par Weierstrafi et pourtant ce sont les methodes qualitatives de Poincare et 
non elle qui nous renscignent sur le comportcmcnt, cventucllement chaotique, du probleme 
des trois corps. 
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Figure 3: (a) La variable complexe z = x + iy et les cercles de convergence et 
les singularites des fonctions (b) fi(z) = 1/(1 — z) et (c) f 2 (z) = 1/(1 + z 2 ). 
Les deux series pour fx et f 2 autour de z = convergent a l'interieur du 
cercle \z\ = 1. 



on peut ecrire en coordonnees polaires z = pe ou 9 est l'angle indique 
dans la figure [3] et p = \z\ = \fzlt est la longueur de la ligne entre l'origine 
(0, 0) et le point z. Une fonction uniforme / est definie en associant a 
chaque nombre complexe z un seul autre nombre complexe f(z). Exemples 
des fonctions uniformes sont des fonctions algebriques comme f(z) = 1+z 
ou (2 + z)/(l + z 2 ) ou des fonctions entieres (c.a.d. analytiques dans le 
plan complexe entier, sauf a Thinni z = oo) comme e z , cos z ou sin z. A cause 
de la relation d'Euler 

e lz = cos z + i sin z (Al) 

les trois dernieres fonctions ne sont pas algebriquement independantes. 
II existent aussi des fonctions multiformes, comme f(z) = y/z, ou a 
chaque valeur de z correspondent plusieures branches de f(z) avec des valeurs 
differentes. 

Comme on le fait pour des fonctions reelles, on peut introduire la notion 
de la derivee f'(z) = df(z)/dz d'une fonction complexe. On dit que / est 
analytique dans un domaine D du plan complexe si pour tout z } zq e D, 
f'(z ) = lim z ^. Zo (f(z) — f(zo))/(z — zq) existe. Si / est analytique dans 
un domaine autour d'un point Zq, il existe un nombre maximal p > et 
des constantes f n tels que f(z) = Y^^Lo fn(z ~ z o) n s'ecrit sous forme d'une 
serie uniformemen^ convergente pour \z — Zq\ < p. p s'appelle rayon de 

33 La definition precise de la notion de convergence uniforme d'une serie comme utilisee 
par Weierstrafi et la difference avec la convergence ponctuelle se trouve dans tous les textes 
d'analyse. 



17 



convergence. De plus, sur le cercle \z — zq\ = p il existe une singularite 
ou / n'est plus analytique. Un resultat profond de la theorie des fonctions 
complexes est qu'on peut caracteriser une fonction analytique par le lieu et 
la nature de ses singularites, voir [Knopp 1971]. 

Nous ne donnons ici que deux exemples illustratifs. Considerons les fonc- 
tions (qui fournissent des exemples des fonctions meromorphes) 

Hz) = ; h(z) = -Lj (A2) 

Si \z\ < 1, on peut representer ces deux fonctions a l'aide des series conver- 
gentes autour du point z = 

oo oo 

/!(*) = £*" ; M*) = EC- 1 )"* 2 * (A3) 

A l'interieur du cercle \z\ = 1 indique dans la figure |3|, les series conver- 
gent uniformement. Le domaine de convergence, qui est l'ensemble des z 
complexes tels que \z\ < 1, est limite par la presence d'une ou plusieures sin- 
gularites sur le cercle \z\ = 1. Pour la fonction fi, cette singularite est visible 
comme une singularite sur l'axe reel en z = 1. En revanche, les singularites 
de la fonction f 2 sont en z = ±i, hors de l'axe reel. Meme si f%{z) n'a aucune 
singularite apparente pour toute valeur reelle de z, ses singularites complexes 
empechent la serie associee de converger pour \z\ > 1. 



Annexe B. Sur les collisions binaires 

La notion de collision binaire dans le probleme des deux corps est presentee de 
maniere elementaire. De part la conservation du vecteur du moment cinetique 
L, le mouvement des deux corps est situe est dans le plan perpendiculaire 
a L. Les positions dans ce plan seront definies par les coordonnees (x, y) et 
la distance entre les deux corps par le module r = \z\ du nombre complexe 
z = x + iy. Apres un changement de coordonnees, l'equation du mouvement 
s'ecrit ^ 2 

dt 2 r 3 ^ 
La propriete de conservation de I'energie s'etablit ainsi : pour z et son 
conjugue complexe zona d 2 z/dt 2 = —z/r 3 et d 2 z/dt 2 = —z/r 3 . En mul- 
tipliant la premiere par dz/dt et la seconde par dz/dt, la somme des deux 
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equations donne 



dz d 2 z dz d 2 z 

did^ + dtdt 2 



1 / dz dz 
■— \z— + z— 



dt 



dt 



(B2) 



Ceci s'ecrit aussi sous la forme 

d f dz dz 
dt \ dt dt 



dt 



(B3) 



Par integration, on trouve 



dz 



dt 



2 (~ + h 



(B4) 



ou la constante h est l'energie totale du mouvement. 

L'equation (|H]) se deduit par changement de variables : Dans la suite, 
nous ecrivons z — dz/dt et z' — dz/ds et de meme pour r et w. En appliquant 



le changement de variables (0) on obtient z 



r z et z 



2ww'. De plus, 



\w\ 2 et la conservation de l'energie ( B4|) s'ecrit dans les nouvelles variables 



sous la forme 



2\w 



l\2 



A l'aide des relations z 
vement ( |B1[ ) s'ecrit sous la forme 



r 3 lrw" — r'z'] 



l + hr (B5) 
et w~ l = w/r, l'equation du mou- 



2w" - (2\w'\ 2 - l) 



1)™ 

r 



(B6) 



L'equation ([10|) se deduit immediatement de ( [Bop et (|B6|) . 

Dans la suite de cette annexe, nous discutons la regularisation du mou- 
vement a travers une collision binaire. On peut toujours admettreQ que le 
mouvement est situe sur l'axe x reel. La distance entre les deux corps est 



decrite par l'equation du mouvement d x/dt 



D'abord on s'interesse 



a trouver la fonction inverse t = t(x) qui satisfait l'equation 



« = (2 + 2 hx) ( * 



(B7) 



34 A l'aide d'un changement de variables qui laisse invariant les equations de mouvement. 
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En integrant directement, on obtient 



( - ( ° = / d *#T2fa < B8) 

ou t est une constante. L'integrale se calcule, pour h ^ 0, a l'aide d'un 
changement de variables x = h' 1 sinh 2 (u/2) et la solution finale s'ecrit (sous 
forme parametrique) 

(2hf /2 (t - t ) = sinh(w) - u , 2hx = cosh(w) - 1 (B9) 

Lors d'une collision, x tend vers 0. Au voisinage de 0, apres l'equation 
(|B"9|), on a x ~ u 2 . De meme, sinhf/u) — u ~ u 3 ~ t — t . On a au voisinage 
d'une collision 

x ~ (t - t ) 2/3 (BfO) 

ou l'on identifie to comme etant l'instant de la collision. 

La serie t — to = u 3 J2^Lq d u u 2u est convergente pour tout \u\ < oo 
et monotone en u pour tout u reel.Q On en deduit done que la fonc- 
tion inverse u = u(t — t ) existe et se developpe formellement en serie 
u — (t — to) 1 ^ 3 J2^=o K (t — to) 2u ^ 3 ou les b u s'obtiennent a partir des a u . En 
injectant ceci dans l'expression de x{u) donnee par l'equation (p9|), on peut 
trouver des coefficients B u tels que 

oo 

x = x(t-t ) = (t- t f 3 £ B v h v (t - t ) 2u/3 (Bf f ) 

Une telle serie a un rayon de convergence non nul si la fonction complexe 
qu'elle represente est derivable. Considerons alors la fonction complexe x = 
x{t) ou r = (2/i)~ 1 (sinh(u) — w) 2//3 . Cette fonction est definie si la fonction 
t = t(u) a une fonction inverse u = u{t). Pour ceci il faut que la derivee 
dr/du ne 'sannule pas. Or, dr/du = si u = u k = 2nik ou k = 0, ±1, ±2, . . . 
est un nombre entier arbitraire. Dans la suite, nous ecrivons = r(uk) = 
(smh(u±k) — u±k)/(2h). Done la derivee 

dx dxdu da; /dr\ 1 3 sinh(M) [sinh(tt) — u] l l 3 (B12) 
dr du dr du\du) 2 cosh(w) — 1 



.35 



Ici, a v = {2h)-^ 2 /{{2v + i)\). 
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existe si u ^ Uk- Autour de uq = 0, elle existe aussi, parce que 



^.ta^y^r 3 (B13) 

n ^o dr u-o 2 V 6 / V 2 / 

est une constante finie. Comme u = correspond a r = et done a t = to, la 
fonction x(t) est derivable au voisinage de la collision binaire et n'a aucune 
singularity dans le domaine |r| < |r±i|. La serie formelle (Bll) a done un 



rayon de convergence p > |r±i| > et represente une fonction analytique. 

Si u = Uk, on a x = a l'instant t donne par r = 7).. Par consequent, les 
Uk caracterisent une suite de (pseudo) collisions. Le passage de t a r nous a 
permis de regulariser les equations afin de pouvoir continuer analytiquement 
le mouvement au-dela la collision en t = to. La non-analyticite apparente a 
T = T k pour fc ^ indique simplement que le choix de r n'est pas encore 
optimal et devrait etre remplace par la variable s definie par Sundman.0 

Notons qu'on a aussi une analyticite de x en fonction de h. 



Annexe C. Un exemple de serie de Lindstedt 

Nous exposons brievement la methode perturbative de Lindstedt, en suivant 
[Rand et Armbruster 1987]. A titre d'exemple, et dans un souci de sim- 
plicite technique, nous considerons l'oscillateur de van der Pol, defini par 
l'equation differentielle non lineaire 

ou x = x(t) est la fonction a determiner en fonction du parametre Pour 
simplifier l'expose, nous poserons comme conditions initiales x(0) = So et 
dx/dt(0) = 0, ou S est une constante. On sait que pour toute valeur de e, il 
existe une solution periodique, appelee cycle limite et que x{t) va toujours 
evoluer vers le cycle limite lorsque t tend vers l'infini. Comme seul le cas 

36 Pour le rapport avec les variables utilisees par Sundman, on a au voisinage d'une 
collision x{t) ~ (t — to) 2 ^ 3 et u ~ (t — to) 1 / 3 , done du = — < )~ 2 / 3 d£. On retrouve ainsi 
une propriete de la regularisation de Sundman dans (^|). 

37 Si e = 0, on retrouve l'equation d'un oscillateur harmonique. 
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e = est facilement soluble, on se propose de resoudre le cas general sous 
forme d'une serie perturbative 



x(t) = x (t) + exi(t) + e 2 x 2 {t) + 



(C2) 



et d'utiliser un changement de variable temporelle r — ut propose par Lind- 
stedt ou 

u = l + k 1 e + k 2 e 2 + ■■■ (C3) 

En ecrivant f{r) pour la derivee par rapport a r de la fonction /(r), on 
obtient le systeme d'equations, jusqu'aux termes du second ordre en e inclus 



Xq + Xq 

X\ + X\ 
II I 

x 2 + x 2 





x \l 
x x \l 



Xp 



2/ciXq 



(C4) 



2x x 'xi — 2kiX\" — (2k 2 + kl)x " + ki(l — xl)x ' 



qui peut etre resolu par recurrence. Les conditions initiates S = X + eXi + 
e 2 X 2 + • ■ ■ se developpent de la meme maniere et on tient comme conditions 
initiales a l'instant t — 



Xi(0)=Xi , Xi'(0) = ; z = 0,l,2,... 

On deduit des equations (|C4|) , (|C5| ) que xq(t) = X cosr et 
rl „ 1 



:X n — Xn 



4" 



sin r + 2kiX cos r + —Xq sin(3r) 



(C5) 



(C6) 



La resolution de (|C6|) peut conduire a l'apparition de termes dits seculaires 
de la forme rsinr ou rcosr. De tels termes croissent indefmiment avec le 
temps et rendent les developpements non-convergents. Pour les eliminer, il 
faut que les coefficients de sinr et cost dans l'equation ( JCJtjj ) s'annulent, d'ou 
X — 2 et k\ — 0. Par recurrence, on obtient 



x 2 + x 2 



4k 2 + 



xi(t) 

1 



1 3 

■- sin(3r) H — cos(3r) + X\ cost 
4 4 



(C7) 



3 5 
cos r + 2X\ sin r — - cos(3r) + 3Xi sin(3r) + - cos(5r) 



L'elimination des termes seculaires amene a poseiQ X\ = 0, k 2 = 
Notons encore que l'amplitude et la frequence ui de la solution periodique sont 



16" 



:',s 



Le fait que ^ demontre que ui ^ 1 est necessaire. 
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completement fixees par la non linearite de l'equation (|C1] ) et ne dependent 
pas des conditions initiales. A partir des expressions explicites on voit aussi 
qu'on obtient le cycle limite sous forme d'une serie trigonometrique dont 
nous avons trouve les premiers termes 

e e 2 
x(t) = 2 cos{ut) + - [3 cos(3a;t) - sin(3u;t)] + ■•• , u = 1 - — H (C8) 
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